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ABSTRACT 


I : FORTRAN computer program RAWINPROC accepts output from NASA 

i. Wallops computer program METPASSl, and produces input for NASA com- 

puter program 3.0.0700 (ECC-PRD). The three parts together form a 

‘\ 

software system for the completely automatic reduction of standard 
RAWINSONDE sounding data. RAWINPROC pre-edits the 0.1-second data, 
including time-of-day, azimuth, elevation, and sonde-modulated tone 
frequency, condenses the data according to successive dwells of the 
tone frequency, decommutates the condensed data into the proper 
channels (temperature, relative humidity, high and low references), 
determines the running baroswitch contact nimber and computes the 
associated pressure altitudes, and interpolates the data appropri- 
ate for input to ECC-PRD. 



INTRODUCTION 


The University of Utah, under sponsorship of the NASA Wallops 
Flight Center, has developed software for the automatic digital 
processing of data transmitted from the standard RAWINSONDE meteo- 
rological sounder. The following describes the software, RAWINPROC 
(Appendix A), which interfaces with Wallops routines NETPASSl which 
reads the magnetic tape from the field to produce the input file of 
raw data, and ECC-PRD which processes the output file of pressure 
(altitude), reference frequency, temperature, and relative humidi- 
ty, all tabulated uniformly at one-minute intervals (Appendix B) . 

After a brief discussion of the purpose and approach, a sub- 
routine-oriented description is presented which closely relates to 
the annotated code (Appendix A). The input deck is described in 
detail under MAIN, and the control card deck is described in Appen- 
dix C. Flow diagrams, narrative description, the CALL list, and a 
complete glossary of variables is included for each subroutine. 
The variable list includes cross-references, descriptions, units, 
constant values, range, limits, and effects where appropriate. 
Files used by RAWINPROC are identified in Appendix D. 

Purpose of RAWINPROC 

The principal purpose of RAWINPROC is to supplant the routine 
manual processing of RAWINSONDE data. Given at one-minute inter- 
vals from a RAWINSONDE sounding, the time-of-day, the reference 
tone frequency, temperature and relative humidity ordinates [Ref. 
1] from the radiosonde, and the pressure (altitude) and azimuth and 


elevation eaglet of the sounding balloon, existing NASA Computer 
Program 3.0.0700 (ECC-PRD) computes and presents the corresponding 
meteorological data in user-ready form. RAWINPROC provides these 
inputs for ECC-PRD without the usual manual processes of reading 
and interpreting the AN/TNQ-5 pen recorder chart. 

Input data for RAWINPROC is provided by NASA Computer Program 
METPASSl which converts field-recorded sonde and tracking data to a 
convenient format. The field-recorded data includes the time-of- 
day, sonde- transmit ted signal (tone frquency) , and tracking angles 
(azimuth and elevation), all sampled at 0.1-second intervals. The 
received signal from the sonde is an audio tone whose frquency (5- 
203 Hz) is determined by the magnitude of the quantity (channel) 
being measured. A baroswitch relects (commutates) the channels in 
sequence (see below) as the balloon rises. Additional input con- 
stants (calibration data, launch time, etc.) are provided in the 
input card deck. 

Radiosonde Commutator 

The four channels transmitted are temperature, relative hu- 
midity, reference, and high reference [Ref. 1]. The reference 
frequency corresponds to a fixed 95 percent of full-scale output 
frequency of the sonde. The frequencies of the temperature and 
relative hmidity channels in ratio with the reference frequency 
provide at the receiver the fraction of 93 percent full-scale read- 
ing of the temperature and relative humidity sensors, independent 
of electronic gain changes during flight. Successive switching 


from channel to channel ia performed by a preasure-actuated commu- 
tator (baroawitch) [Ref. 1, p. B2-*19] . The leading edge of auccea- 
aive reference, high reference, and relative huaidity contacta of 
the baroawitch correapond to calibrated preaaure altitudea. Tem- 
perature ia transmitted betrnen each of theae contacta. At high 
reference contacta reference frequency is increased a few hertz to 
distinguish it from other reference contacts so as to eliminate 
ambiguity in associating switch points in the received signal with 
their baroawitch contact nunbers. 

The radiosonde commutator bar is represented in Fig. 1. The 


Increasing altitude. 



Contact number 

Switch points 

Baroswltch 
contact pads 


Channel 


Fig. 1. The baroawitch wiper arm traverses the commutator bar as 
the atmospheric pressure decreases during balloon ascent. 
Electric contact with each channel, reference (R), high 
reference (P) , relative hmidity (H) , is made as the wiper 
passes over each commutator pad. Temperature (T) is trans- 
mitted between pads. "Switch points" in the data occur at 
the leading edges of the pads. 
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entire stenderd baroswitch sequence is displayed in Table 1. As 
pressure decreases, contact msnber increases. Baroswitch output 
dwells on relative hunidity, reference, and high reference signals, 
according to successive contacts as tabulated. Temperature data 
are transmitted during intervals between contacts. The leading 
edge of each of the contacts represents a pressure, calibrated for 
each radiosonde. The contacts are traversed at a rate depending on 
the balloon rise rate. Spacing between contacts approximately 
equals contact width, so baroswitch temperature dwells are of 
length comparable to contact (humidity, reference, and high refer- 
ence) dwells. 

The baroswitch output, then, alternates between temperature 
and either relative humidity, reference, or high reference, depend- 
ing on contact number. Contacts below number 135 transmit relative 
humidity (H), except every fifth contact. The fifth contacts 
transmit reference frequency (R), except that every third reference 
beginning with contact number 30 is high reference (P). 

Beginning with contact number 135, no humidity is transmit- 
ted. Each contact, H5 to 179, transmits reference frequency, with 
every fifth one a high reference. The pattern is recognizable in 
the received data so that pressure (sonde altitude) can be assigned 
to the common time base of the data channels. 

Decommutation Approach 

The approach taken in DECM is first to track the input fre- 
quency function from the sonde and to detect commutator switching 
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Table 1. BAROSWITCH SEQUENCE 


Contact Nimber and Channel Tiannitted 


1 2 3 4 5 6 7 
H H H H H H 

33333333 
0 1 2 3 4 5 6 7 
H H H H „ H H 
P 

66666666 
0 1 2 3 4 5 6 7 
H H H H . H H 
P 


99999999 
0 1 2 3 4 5 6 7 
H H H H „ H H 


11111111 

22222222 

0 1 2 3 4 5 6 7 

H H H H „ H H 
R 


11111111 

5 5 5 5 5 5 5 5 

0 1 2 3 4 5 6 7 

„RRRR„RR 
P P 


11111 
8 9 0 1 2 3 4 
H H H H H H 

3 3 4 4 4 4 4 
8 9 0 1 2 3 4 
H H H H H H 


6 6 7 7 7 7 7 
8 9 0 1 2 3 4 
H H H H H H 


11111 
9 9 0 0 0 0 0 
8 9 0 1 2 3 4 
H H H H H H 


1111111 
2 2 3 3 3 3 3 
8 9 0 1 2 3 4 
H H H H H H 


1111111 
5 5 6 6 6 6 6 
8 9 0 1 2 3 4 
R R R R R R 


111112 2 

5 6 7 8 9 0 1 

, H H H H „ H 

R R 

4 4 4 4 4 5 5 

5 6 7 8 9 0 1 

H H H H _ H 

P 

7 7 7 7 7 8 8 

5 6 7 8 9 0 1 

H H H H „ H 

P ^ 

1111111 
0 0 0 0 0 1 1 

5 6 7 8 9 0 1 

H H H H „ H 

P 

1111111 
3 3 3 3 3 4 4 

5 6 7 8 9 0 1 

„ R R R R _ R 

P P 

1111111 

6 6 6 6 6 7 7 

5 6 7 8 9 0 1 

„ R R R R „ R 

P P 


22222222 

23456789 

5 5 5 5 5 5 5 5 
23456789 


88888888 
23456789 
HHH HHHH 


11111111 
11111111 
23456789 
HHH HHHH 


11111111 

44444444 

23456789 

RRRpRRRR 

11111111 
77777777 
23456789 
RRR RRRR 


Read: contact nmber 13:H, 165:P, etc., represented in the 

1 1 

table as 3, 6, , etc. 

H 5 
P 

H: relative humidity 

R: reference 

P: high reference 

(T, temperature is transmitted between each contact pad) 


events. The dete ere inmediatelyi Chen, condensed to the mean 
frequency and midpoint time (one datum) for each signal dtiell. The 
switching time at the leading edge of each dwell, and the length of 
dwell are also stored. Tracking frequency gates as well as tine 
gates are used to identify the channels. "Guard sones", first- 
order extrapolation, and other added logic are used to reduce sus- 
ceptibility to channel crossovers, sudden changes of signal in a 
given channel, and to variabilities in sonde design and perfor- 
mance. 

Third, the program assigns contact nunbers, and therefore 
pressures, to the appropriate signal switch times according to the 
baroswitch sequence. 

Finally, the desired output table at one-minute intervals is 
constructed by interpolation from the asynchronous decommutated 
time functions. 



SUBROUTINE DIAGRAM OF RAWINFROG 









RAWINPROC(MAIN) 


D«»cription. MAIW 

HAIN includes the segments of code described ss Initialiser, 
Advancer (S.2S4), Condenser (S.265, Oecommuta^ci and Baroswitch 
Tracker (S.410), and Terminator (S.81) in the block diagram (Fig. 
2). The present discussion describes each of these segments and 
presents a detailed description of the input card deck for 
RAWINPROC and a glossary of all the variables in MAIN. Separate 
discussions below describe each of the subroutines ADVANCE, ANCLE, 
TRACK, SEARC', DECOM, and INTERP, includir« the flow through th .> 
respective CALL lists. The discussion closely follows the FORTRAN 
program list (Appendix A). Throughout this document zeroes which 
might be mistaken as letter *'0" are given the slash, 0 . Also, the 
term Humidity is used in place of the longer Relative Humidity. 

Initializer 

Input parameters are read from the punched card input deck 
(see below) and printed. 

Launch time (hours, minutes, seconds, GMT) is read in, to- 
gether with optional limiting elapsed time of processing (TPROC), 
and optional beginning time of processing (TSTART). When TPROC is 
unpunched in the input deck, processing terminates on other 
criteria. When TSTART is unpunched, processing starts two minutes 
before balloon release (or at the beginning of recorded data, 
whichever is later). 









Elonents of the COND ( , ) array are printed by MAIM, if 

desired for diagnostic purposes, along with other information 
indicative of processing progress. Page headings are printed every 
LINENO lines. Additional optional printout is available for 
diagnostic purposes according to input values assigned to TbST(I), 
where 1 * 1-3 is used in ADVANC, SEARCH, and TRACK; 1 * 4, 5 in 
SEARCH; 1 ■ 6 in INTER? ; and 1 ■ 7-9 in MAIN, 

Condenser variables FSUM, MSUM, LOS, and JK are initialized, 
and constants HGATE - 1.0 (Hz), SIGMIN « 5.0, SIGMAX » 205.0, and 
IN are assigned or computed. The COND/lCONT) array, the raw data 
counter JJ, and the one-minute table VL are all cleared. 

OUTPUT variables TNOH, LIST, and ISTOP are initialized, 
DLIST, TGHDAQ, and PCAL are read in, and LCNTK defined. Obvious 
errors in PCAL are automatically corrected and identified. Surface 
readings FP0, FTEMP0, FRH0, and FR0 are read in, converted where 
necessary, and stored as initial values of V2. The decimal contact 
number (AICR0) equivalent to surface pressure is computed and 
stored in integer form (ICR0) for use by DECOM. 

Manual burst input CBRST is read in, and TBRST is initialized 
for DECOM. 

The initial frequency TF is computed from the surface ordi- 
nate readings according to 

f » e * 2. * FR0/95.O 

where f represents the frequency corresponding to ordinate 6, and 
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FR0 is Che ordinate reading of low reference signal (at the time of 
launch) when the AN/TMQ-S pen recorder has been adjusted to 30 
ordinates for 60 Hz input. Assuming recorder linearity the ratio 
of frequency to ordinates is the same at reference frequency fj^ as 
at 60 HZ) so 


—7 * ~ . f « 2*FR0 
FRff 30 ’ R ^ 

Since, after obtaining FR0, the recorder is continuously adjusted 
to maintain 95 ordinates at reference frequency, then 

f h 

^ 95 

Thus, for the initial reference frequency fj^ ■ 2*FR0, the first 
expression above converts initial ordinate values FTEMPO and FRH0 
to their corresponding sonde frequencies TF and HF. The resulting 
initial frequency f (i.e., TF) is used by OECOM to initially 
position frequency gates. (HF has proved not useful in DECOM.) 

Subroutine ADVANC is initialized by reading Che raw data from 
the beginning of the file (FILE 01) , Appendix D) until the forward- 
most point, TIME(IO), of a ten-point sample is past TSTART. To 
protect against time word errors ("spikes" or constants), the fol- 
lowing conditions were imposed before accepting the starting point: 

-1.0 < TIME(IO) - T1ME(6) - 0.4 < 1.0 seconds 
TIME(IO) - T1ME(6) >0.2 seconds 
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or, equivalently: 


ITIHE(IO) - TIME(6) - 0.8| < 0.6 


Advancer 

The raw data file, TIME, FREQ, AZ, EL, at 0.1**second recorded 
data rate, is processed by Condenser ten points at a time. The 
sample of ten points, however, is advanced only five points at a 
time. Condenser therefore searches for and tracks signal and de- 
tects switch tines between signal dwells by examining in sequence 
half-overlapping 1.0-second samples of raw data. At each return 
for more data. Advancer moves the 1.0-second ten-point sample 
(TIME, FREQ, AZ, EL) ahead one-half second. Advancer also counts 
the mmber (JJ) of raw points read in, for diagnostic purposes, and 
causes termination if TSTOP is passed. 

In addition, at each one minute after launch, except before 
AN/GMD acquisition time TGMDAQ, Advancer sends the ten-point sample 
to subroutine ANGLE to compute the output values of AN/GMD angles 
AZ and EL. The latter t%» quantities are stored for subroutine 
INTER? in VL(2, ) and VL(3 ), along with the associated elapsed 
minutes from launch in VL(1, ). Though this processing of angles 

is a condensing function, it is nevertheless more conveniently 
located in the code of the Advancer. 

Condenser 

Condenser (subroutine TRACK) determines whether the ten-point 
sample lies sufficiently in the signal-tracking gate. If so TRACK 
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adjusts the gate, accumulates data toward the mean frequency of the 
current d%iell, and returns for the next ten-point sample, repeating 
this process until the signal does not lie sufficiently in the 
tracking gate. Vfhen the signal changes abruptly, subroutine TRACK 
calls subroutine SEARCH to reposition the gate on the signal. 
SEARCH returns via TRACK to MAIN (advancer), (JR ■ JKMEM) for new 
data until it finds signal. When signal is found and the tracking 
gate (SIGLEV) is repositioned, a new condensed point (JR ■ JR 1) 
is stored in COND( , JR), representing the preceding signal 
dwell. Upon returning to MAIN, control proceeds to the Decommuta- 
tor and Baroswitch Tracker (subroutine DECOM) . 

It is noted that the production mode (ITYPE ■ ”P”, using no 
File 03) is the principal and normal mode of operation of 
RAWINPROC. Other modes were used during program development for 
economical reasons as follows: 

ITYPE ■"M" includes the wtiting of File 03, Condenser 
output file 

ITYPE »"C" reads File 03 (COND, VL) , skips Advancer, 
Condenser 

Decommutator and Baroswitch Tracker 

Decommutation of a condensed point COND( , JR), i.e., deter- 
mination of its channel ICONDd, JR) %«hether Temperature, Refer- 
ence, High Reference, or Humidity, is performed by subroutine 
DECOM. The determination of its contact manber, IC0ND(2, JR), when 
it is a Reference or Humidity contact switch point, is also per- 
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foraed by DECOM. Finally, DECOH senses and computes the tine of 
balloon burst, TBRST. 

Terminator 

RAWINPROC discontinues processing raw data upon one of the 
following conditions: 

1. TBRST encountered (MAIN S.410, S.84) 

ISTOP - 10 (INTERP S.2003, S.47) 

2. End of raw data. File 01 EOF (MAIN S.40, S.82) 

ISTOP « 6 

3. Reached TSTOP (MAIN S.81, ADVANC S.3) 

ISTOP ■ 7 

4. COND array overflow, JK > 1000 (MAIN S.85, TRACK S.66) 
ISTOP ■ 8 

5. Excessive loss of signal, LOS > 100 (MAIN S.83, TRACK 
S.83, SEARCH S.IO), ISTOP • 5 

Upon arriving at the terminal exit, S.90, terminator decodes 
ICONDd, ) and IC0ND(2, ) (eliminating diagnostic information 
coded by DECOM). An auxiliary listing of the asynchronous output 
of DECOM is printed if requested by input quantity TEST(7). Then 
subroutine INTERP is called to compute and interpolate VL(I, L) , 1 
* 4-7, pressure (mb), reference (He), temperature (ordinates), and 
relative hianidity (ordinates), at one-minute intervals correspond- 
ing to VL(I, L) , I ■ 1-3, elapsed time and tracking angles 
(azimuth, elevation) tabulated by ADVANC. 

Terminator prints the one-minute table, VL(I, L), L ■ 1-LIST, 
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and terminating comments according to ISTOP. Apogee rows of VL 
having aero pressure are eliminated, as the input "deck" (File 02) 
is prepared for ECOPRO. This terminates activity No* 2 of the 
RAWINPROC system. 













Flow Diagram, MAIN (continued) 











Input Card Deck 


Card 

Column 

Format 

Variable 

Value 

Unite 

Comments 

1 

1-2 

12 

lOIN 

5 


Input deck file 
number 


3-4 

12 

lOOT 

6 


Print file number 


7 

11 

ITYPE 

P 


Production run 

2 

1-2 

12 

11 

(a) 

hour 

Launch time, GMT 


3-5 

13 

12 

(a) 

minute 

Launch time, GKT 


7-10 

F4.1 

TS3 

(a) 

second 

Launch time, GMT 


11-20 

Flo. 2 

TPROC 


seconds 

(d) 


21-30 

F10.2 

TSTART 


seconds 

(d) 

3 

1-5 

F5.0 

TEST(l) 

0. 


(b) 


6-10 

F5.0 

TEST(2) 

0. 


(b) 


11-15 

F5.0 

TEST(3) 

0. 


(b) 


16-20 

F5.0 

TEST(4) 

0. 


(b) 


21-25 

F5.0 

TEST(5) 

0. 


(b) 


26-30 

F5.0 

TEST(6) 

0. 


(b) 


31-35 

F5.0 

TEST(7) 

0. 


(b) 


36-40 

F5.0 

TEST(8) 

0. 


(b) 


41-45 

F5.0 

TEST(9) 

0. 


(b) 


46-50 

F5.0 

TEST(IO) 

0. 


(b) 

4 

1-10 

FlO.l 

DLIST 

60. 

seconds 

(c) 60. required 


11-20 

FlO.l 

TGMDAQ 

20. 

seconds 

Recom. 20. 

5 

1-10 

FlO.l 

FP0 

(a) 

mb 

"SFC PRESSURE” 


11-20 

FlO.l 

FTEMP0 

(a) 

ordinates 

"TEMPERATURE" 


21-30 

FlO.l 

FRH0 

(a) 

ordinates 

"RH" 






6 

7-29 

30 


31 


32 


Column 

Format Variable Value Units 

Comments 

31-40 

FlO.l 

FR0 

(a) ordinates 

"UNADJUSTED 

ORDINATE" 

1-10 

F10.2 

CBRST 

(a) 

(f) 


8F10. 

1 PCAL(i) 

, (e) mb 




i - li 

180 


2-17 

45 (DATE ATOOAYAyyddd (yy ■ year, ddd ■ day of year)] 

21 

11 

(JOPT) 

2 

<g) 

25 

11 

(KOP) 

0 

(g) 

28-29 

12 

(IN) 

5 

(8) 

30-31 

12 

(10) 

6 

(g) 

32-33 

12 

(IT) 

9 

(g) 

36-37 

12 

(WOPT) 

(blank) 

(g) 

75-80 

A6 



l.D. , ascent 
number, e.g., 
AS607A 

2-24 

A6 (THISAISAAAECCAAARAWINSONDE) 

(g) 

31 

11 

(lAZ) 

1 or 0 

(h), (g) 

33-39 

F7.5 

(DOBSON) 0.0 

(g) 

75-80 

A6 



l.D. , ascent 
mrnber, e.g., 
AS607A 

2-11 

2A5 

(ISTTI, 

1STT2) [WALLOPSAIS] 

(g) 

13-18 

17 

(LDATE) 

(nmddyyl 

(g) 


(mm ■ 

month, dd 

■ day of month, yy 

■ year) 

20-25 

17 

(UNCH) 

(hhmm72] 

(g) 


(hh ■ 

hours, mm 

■ minutes, (GMT), 72 ■ WALLOPS) 

26-31 

F6.1 

(HTMSL) 

4.0 meters 

(g) 

34-37 

F6.1 

(CALTP) 

30.0 

(g) 


- 19 - 
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Input C«rd D«ck (continued) 


Colunn 

Fonut 

Variable 

Value Units 

CosBMnts 

39-43 

F6.1 

(RECTP) 

(a) 

Ordinate value 


corresponding to 
30*C (Tenporsture 
Calibration) 


45-49 

F6.1 

(CALRH) 

(a) 


Ordinate value 
corresponding to 
-40* and 46 ord. 
(im calibration) 

52-55 

F6.1 

(RECRH) 

46.0 


(g) 

56-57 

12 

( ICBRN) 

02 


(g) 

58-62 

F5.1 

(SORT) 

(a) 

*C 

Surface 

temperature 

63-67 

F5.1 

(SI!<H) 

(a) 

X 

Surface RH 

69-74 

F6.1 

(PCAL) 

(a) 

mb 

Surface pressure 

75-80 

A6 

(ID) 



l.D. , ascent 
mnber, e.g., 
AS607A 

2-11 

A6, A4 (NLBL(2) 

,( 3)1 WALLOPS AIS) 

(g) 

13-18 

A6 

(nwi ■ 

(DLBL(l) 
month, dd 

(mmddyyl 

• day of month, yy 

(g) Launch data 
■ year) 

20-69 





No osonesonde 
calibration need- 
ed 

75-80 

A6 




l.D. , ascent 
ntanber, e.g., 
AS607A 

2-7 

17 

(mm ■ 

( ILDTE) 
month, dd 

(naddyyl 

■ day of month, yy 

(g) Launch data 
■ year) 

9-12 

17 

(hh - 

(ILTME) 
hour, sn ■ 

thhmm72l 

minute. 

(g) Launch time 
ss ■ atation number ■ 72) 

15 


(minus s 

ign) 1-) 


(g) 

16-19 

15 


19999] 


(g) 
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Input Crd Ptck 


Column 

Format 

Variable 

Value 

Units 

CoMents 

21-26 

P7.1 

(HCMDT) 

4.0 

meters 

(g) Ceopotential 
height of station 

29-31 

F5.0 

(VSFC) 


meters/s 

(g) Surface wind 
speed 

34-36 

F5.0 

DSFC 


degrees 

(g) Surface wind 
direction 

39-74 






75-80 

A6 




I.D., ascent 
number, e.g. , 
AS607A 


Comments: 

(s) Derived from field launch records. 

(b) Other values used only for internal test purposes. 

(c) Obsolete input, vslue fixed. 

(d) Normally unpunched (sero). 

(e) Card deck punched from the baroswitch pressure calibra- 
tion data (chart or punched paper tape) provided with 
each sonde. 


(f) The terminating baroswitch contsct ninber, CBRST, to the 
nearest one-hundredth the distsnce between switch points 
(contact leading edge) , is provided by standard manual 
procedure (Ref. 1, pp. B4-B12, BIT, BS-B8). Termination 
by CBRST is allowed only after 50 minutes (30QC seconds) 
of flight. Such termination also may be useful for s 
variety of other reasons, such as battery or other in- 
flight failure. 

(g) For encoding format of cards 30-34, see ECOPRO documen* 
tation, NASA Computer Program 3.0.0700, NASA Uallops 
Cofliputer Program Abstracts, Vol. XXVll (sonde 10, rsdio- 
sonde and osonesonde calibration, and flight-end cards). 

(h) Input lAZ is determined by which tracking system is 
used. lAZ ■ 0 if asimuth zero is north, ■ 1 if south. 
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Lilt of Vari«bl«i, MAIW 


AlCno Real fora of the initial value of XCR0. 

AZ(IO) Ten-point aanple of O.l-second raw data* aaiieuth 

angle. 

AZK Aainuth, uaed in writing VL(2, ) into ECC-PRO 

input file (File 06). 

CBRST Contact number, to neareit 0.01, at balloon bunt, 

input card No. 6 (lee Input Card Deck, above, and 
•ubroutine DECON). 

CNVOF Conversion factor, ordinates to frequency (Hs). 

C0ND(3, 1000) Condensed point array, real: 

CONDd, ) ■ Time (seconds) from launch 

C0ND(2, ) ■ Duration (seconds) of the signal dwell 

CONOO, ) ■ Mean signal frequency (herta) over the dwell. 

DIFFl, Differences used in detecting and correcting gross 

D1FF2, errors in baroswitch calibration table, PCAL. 

DIFFAV, 

OlFFHl, 

DIFFLO 

DLIST ■ 60 seconds, data rate of ECC-PRO input cards, 

used in computing output array VL(7, LIST), input 
card No. 4. 

Dummy variables, used in reading raw data file, 
(File 01). 


DM1 

DM2 

DUM(IS) 

EL(10) 


Ten-point sample of 0.1-second raw data, elevation 
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ELPT 

ELV 

IH 

FR0 


FREQ( 10) 
FRH0 

FSUM 

FTEMP0 

HF 

HGATE 

I 

II 
12 
IBC 


angle. 

Tine (seconds) from launch of condensed points, 
used in auxiliary listing of OECOM output. 

Elevation angle, used in tnriting VL(3, ) into ECO 

PRD input file (File 06). 

Surface atmospheric pressure (mb) at balloon re- 
lease, input card No. 5. 

Stripchart reading (ordinates) of Reference channel 
when recorder gain is adjusted to 30 ordinates for 
60 Hz test input. Input card No. 5. Used in com- 
puting CNVOF. 

Sample of ten 0.1-second raw data points (Hz). 
Stripchart reading (ordinates) of Humidity channel 
at balloon release, input card No. 5. 

Running sun (Hz) of signal means in signal tracking 
gates for computation of the signal dwell mean 
C0ND(3, ) - FSUM/NSUM. 

Stripchart reading (ordinates) of Temperature chan- 
nel at balloon release, input card No. 5. 

Signal frequency equivalent to FRH0 (not used). 
Half-width (Hz) of signal tracking gate, constant 
1.0 Hz. 

DO-loop index. 

Hours integer of launch GMT, input card No. 2 
Minutes integer of launch GMT, input card Mo. 2. 
Intermediate variable used in decoding 1C0ND(2, ) 
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after termination (MAIN S.90) of DECOM. 

IBT ECC-PRD ozone input datum (OZN4) , DNA, IBT > 0. 

IC0ND(2, 1000) Condensed point array, integer: 

ICONOd, ) ■ channel (1 ■ temperature, 2 ■ reference, 3 ■ high 
reference, 4 ■ relative humidity) 

IC0ND(2, ) <B baroswitch contact number 

ICR0 Effective baroswitch contact number at launch (to 

nearest ttio decimal places multiplied by 100) com- 
puted from surface pressure and PCAL table. 

ID ECC-PRD input, quantity (ID), identification, input 

card No. 32. 

IN The number of subdivisions of the sonde frequency 

range used in searching for signal. 

INDAX Working channel number, ICONDd, ) before decod- 

ing, used in auxiliary printout of DECOM output. 

10 ECC-PRD input file number. File 02, 10 ■ 02. 

100 ECC-PRD ozone input datun (0ZN2) , DNA, 100 ■ 0. 

103 ECC-PRD ozone input datun (DJP) , DNA, 103 “ 0. 

IOC ECC-PRD ozone input datum (0ZN3), DNA, IOC * 0. 

lOlN Input card deck file number (File 05), lOlN * 5. 

lOUT Output print file number (File 06), lOUT • 6. 

IPRINT Count of printed input raw data samples, TIME(J) , 

FREQ(J), J » 1, 10, in auxiliary printout; IPRINT ■ 

1 - 10 . 

IRH ECC-PRD Humidity input datum (OF), IRH ■ VL(76, ) 

* 10 . 0 . 
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1ST 


Additive term %ihich increases with altitude the 


I STOP 

ISTTl, 1STT2 
IT 

ITl, 1T2 

IT3 

ITP 

ITY 

ITYPE 

IX, lY 

J 

JC 

JJ 


percentage limit within which PCAL values are test- 
ed in sequence by contact. Used in the process of 
detecting gross errors in the PCAL table. 

Variable indicating termination status (see MAIN, 
Description, Terminator) . 

ECC-PRO input (iSTTl, ISTT2) , station name, input 
card No. 32. 

ECC-PRD input datun (TIM), IT - VL(1, )/60.0. 

Time of day, hours, minutes; used in auxiliary 
printout of DECOM output. 

Time from launch, minutes; used in auxiliary print- 
out of DECOM output. 

ECC-PRO Temperature input datum (DT) , ITP • 
VL(6, )*10.0. 

First character, V, of a ”VL" record in COND/VL 
file. 

Input character (input card No. 1). ITYPE ■ P. 
(See MAIN, Description, Condenser.) 

DO-loop indices, used in reading and printout of 
baroswitch pressure calibration, PCAL, table. 
DO-loop index, used in initializing VL, indexing 
raw data, and testing for zero pressures in VL. 
DO-loop index used in auxiliary printout of DECOM 
output. 

Count of raw data records read. 
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JK 

JKMEM 


JKT 

JP 

KNTCT 

LCNTK 

LDATE 

LINE 

LINENO 

LIST 

LL 

LOS 

LTIME 


Count, or current index, of condensed point COND. 
Value of JK upon each entering of the condensing 
process, used to detect %ihether subroutine TRACK 
requires a new raw data point or has concluded a 
new dwell. 

Minute from launch, integer stored in File 03 with 
associated VL. DNA in production (ITYPE ■ P) runs. 
DO-loop index, used in computing 1CR0. 

Last contact processed by INTERP, used only in 
terminal printout if INTERP encountered LCNTK. 

The highest number contact pressure-calibrated. 
(See PCAL.) 

ECC-PRD input (LDATE) , input card No. 32. 

ECC-PRD input character string, used in reading 
input cards No. 30-34. 

Lines printed per page, used in labeling auxiliary 
printout of DECOM output (File 06). 

Count, or current index, of rows of VL array, i.e., 
manber of one-minute "cards" input to ECC-PRD. 
DO-loop index, used in auxiliary printout of VL 
array and in eliminating zero-pressures at apogee 
in VL. 

Loss of signal count from Condenser (see subroutine 
SEARCH) , used in terminal printout when processing 
terminated due to signal noise (ISTOP • 5). 

ECC-PRD input (LANCH), input card No. 32. 
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NSUM 

PCALdSO) 

FERC 

PR 

RTl 

RT2 

RT3 

RT4 

SIGLEV 

SIGMAX, SIGMIN 

TBRST 

TEST(IO) 

TF 


Quantity initialiaed for Condenaer (aee aubroutine 
TRACK) . 

Calibrated pressure values, in order corresponding 
to baroswitch contact number, input cards No. 7-29. 
Limiting percent change, usdd in automatic checking 
successive PCAL values for gross errors. 

ECC-PRD input (PR), PR « VL(4, ). 

Seconds part of time of day (ITI, 1T2, RTl), used 
in auxiliary printout of DECOM output. 

Time of day (hours), used in auxiliary printout of 
DECOM output. 

Seconds part of elapsed time from launch, used in 
auxiliary printout of DECOM output. 

Same as C0ND(2, ), dwell (seconds), used in auxil- 

iary printout of DECOM output. 

Center frequency (Hz) of signal tracking gate in 
Condenser (see subroutine TRACK). 

Upper and lower limits of sonde frequency range 
(Hz). SIGMIN « 5, SIGMAX- 205. 

Computed time (seconds from launch) of balloon 
burst (see subroutines DECOM, INTERP) . 

Input constants which control diagnostic printout 
(see MAIN, Description, Initializer). 

Initial signal frequency of Temperature channel, 
computed from initial temperature ordinate PTEMP0 
and sonde frequency calibration FR0. Passed in 
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TGMDAQ 


TIME(IO) 


TLANCH 


TNOH 


TPROC 


TS3 


TSTART 


TSTOP 


V2(7) 


VL(7, 150) 


first call to subroutine DECOM. 

Time delay (seconds) after balloon release of ex- 
pected AN/GMD antenna acquisition of the sonde. 
Input card No. 4. 

Time of day (seconds) in the ten-point sample of 
raw data processed by Condenser. 

Time of day (seconds) of balloon release 

(launch). Computed from data on input card No. 2. 
Time (seconds from launch) of the occurrence of 
contact 135, when transmission of humidity data 
ceases. 

Maximum time interval (seconds) of flight data to 
be processed, automatically made large if unpunched 
on input card No. 2. 

Seconds part of time of day at launch, input card 
No. 2. 

Time from launch (seconds) to begin processing 
flight data. Made -120 seconds (tw<} minutes before 
launch) if left unpunched on input card No. 2. 

Time (seconds from launch) at which processing is 
terminated if requested via input TPROC. TSTOP ■ 
TPROC - TSTART. 

Initial values for interpolation in subroutine 
INTER? (see). Computed from surface input quanti- 
ties FP0, FTEMP0, FRM0, FR0. 

One-minute array, output of INTERP, input data for 
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ECC-PRD. 


VL(1* ) ■ Time (seconds from launch) at one-osinute inter* 

vals. 

VL(2| ) • Azimuth (degrees) 

VL(3, ) * Elevation (degrees) 

VL(4, ) ■ Pressure (mb) 

VL(5, ) ■ Reference frequency (Hz) 

VL(6, ) • Temperature (ordinates) 

VL(7, ) • Relative humidity (ordinates) 

XM ECC'PRO ozone input quantity (XM) , DNA, XM • 0.0. 
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SUBROUTINE AOVANC 


Description 

The ten-point sample of 0.1-second raw data, TIME( J) , 
FREQ(J), AZ(J), EL(J), J ■ 1, 10, is advanced five points at each 
call of subroutine ADVANC. A running sin, JJ, of raw data points 
read in is kept by ADVANC for diagnostic purposes only. 

The raw data word TIME is converted from time of day (hours) 
to elapsed time from launch (seconds). The raw data word FREQ is 
converted from period (milliseconds) to frequency (hertz). Those 
with periods outside the range 4.8 to 200 ms (5 to 208 Hz) are made 
zero hertz, to prevent dividing by possible extreme values. If the 
TIME word does not increase by 0.1 ± 0.05 seconds, its value is 
replaced by the preceding value increased by 0.1 second. This 
"corrects" possible "spikes" and other temporary time word errors, 
relying on the fact that the data vere recorded in the field at 
real-time 0.1-second intervals. Gross faultiness in the time word 
is usually detected in MAIN when searching for TSTAPT. 

Vfhen the accepted time reaches TSTOP, ADVANC causes termina- 
tion via RETURNl (MAIN S.82). 

At each minute after field tracker target acquisition time 
TGMDAQ, ADVANC calls subroutine ANGLE to compute from the ten-point 
sample AZ(J), EL(J), J ■ 1, 10, a value each for azimuth and eleva- 
tion to place in the RAWINPROC output table VL(2, ) and VL(3, ), 

respectively. 

If input TEST(l) is greater than zero, and TlME(l) lies be- 
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tween input values TEST(l) end TEST(2)» the angle raw data AZ(J), 
EL(J)( J ■ 1, 10 are printed out* 



Advance later half to earlier half of 
ten-point aample 

TIME(J), FREQ(J), AZ(J), EL(J), J - 1, 10 


Read five points of new raw data into 
later half of the ten-point window. 

a. Convert TIME to seconds from launch. 

b. Confirm/substitute O.l-s time incre- 
ments. 

c. Convert FREQ to hertz. 

d. Substitute zero for FREQ values out 
of range. 


Enter angle data 
into output arra? 


ANGLE I 
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CALL Lilt, ADVANC (Ref.: List of Verieblee, beloti) 


Variable 

Flow 

Cosmwnts 

TIMECIO 

Initial 5 pointa from 
MAIN ( Initial iaer) to 
MAIN (Condenser) 

ADVANC moves this sample 
forward S pointa each call. 

FREQ(10) 

Initial 5 points from 
MAIN (Initialiser) to 
MAIN (Condenser) 

ADVANC moves this sample 
forward S points each call. 

AZ(IO) 

Initial 5 points from 
MAIN (Initialiser to 
MAIN (Condenser) 

ADVANC moves this sample 
forward 5 points each call. 

EL(IO) 

Initial 5 points from 
MAIN (Initialiser to 
MAIN (Condenser) 

ADVANC moves this sample 
forward 5 points each call. 

JJ 

From MAIN to MAIN 

ADVANC increments at each READ. 

TSTOP 

From MAIN 

Used to terminate processing. 

TLANCH 

From Main 

Used in converting to elapsed 
time. 

TGMDAQ 

From Main 

Used to prohibit processing 
meaningless angle data. 

TEST( 1 ) , 
TEST(2) , 
TEST(3) 

From Main 

Control diagnostic printout. 
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Lilt of Variablo, ADVAHC 


AZ(IO) Asimuth angle (degrees), ten-point senple of 0.1- 

second raw data. 

DLIST Time interval (60 seconds) of uniform output table 

VL (see input deck, MAIN). 

DM1, DM2, Dim Dummy variables used in reading raw data file. 


EL(IO) 


FREQ( 10) 


ITYPE 


J, J5, Jl 
JJ 

LIST 


Elevation angle (degrees), ten-point sample of 0.1- 
second raw data. 

Signal frequency (Hz), ten-point sample of 0.1- 
second raw data. 

Program development input to permit rerunning from 
Condenser output file (File 3). Options: 

*'M" writes File 3 (COND/VL array) 
reads File 3 

"P” Production (File 3 omitted). 

Raw data variable indices. 

Total number of raw data points read from File 1. 
Total number of entries stored in one-minute table 


TEST( 10) 
TGMDAQ 


TIME(IO) 


TLANCH 


VL. 

Input variable controlling diagnostic printout. 
Input constant, mnber of initial seconds angle 
data are presumed invalid due to target acquisition 
time required by the AN/GMD balloon tracker. 

Time (seconds from launch) of the ten-point sample 
of 0.1-second raw data. 

Time of day (seconds) of balloon release. 
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TSTOP 


Tim (icconds from launch l«si TSTART) preset by 
card input to stop processing data. 


VL(7| ISO) RAWINPROC output array. The first three colunns 
VL(1, } > (isinute), VL(2, ) - AZ, VL(3, ) • EL, 

are computed and loaded by ADVANC, the remainder by 
INTERP. 
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SUBROUTINE ANCLE 


Detcription 

Subroutine ANCLE edits, condenses, end smooths input AN/GMD 
angle date AZ(1), EL(I), 1 ■ 1, 2, 3, ..., 10 (asinuth, eleva- 
tion) . It computes (when called every DLIST seconds) from the ten 
local consecutive (0.1-second) values, one value assigned at the 
midpoint, AZ(S), EL(S). The computed value is the mean of those 
points lying in the five-degree interval centered on the unit de- 
grees mode of the ten input points. Other points, including ex- 
treme values, are therefore rejected. 

The unit degrees mode is the most populated one-degree inter- 
val over the ten input points. It is found by rounding to units 
place the input values and counting equal rounded values. Vftien the 
distribution is such that more than one unit degree interval has 
the highest population, the one occurring earliest in time within 
the 1. 0-second sample is used. 









CALL List, ANGLE (Ref.: List of Variables, below) 



Flow 

Conments 

AZ(IO) 

From ADVANC 

Ten 0.1-second points, raw 
data, input. 

(AZ(5)1 

To ADVANC 

One 1.0-second condensed point 
output. 

EL( 10) 

From ADVANC 

Ten 0.10-second points, raw 
data, input. 

lEL(5)l 

To ADVANC 

One 1.0-second condensed point 
output. 
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List of Variables, ANGLE 


AZ(IO) 

EL( 10) 

lAZ(IO) 

lEL(lO) 

K 

KMA 

KME 

L 

LL 

NIAZ(IO) 

NIEL(IO) 

NSUMA 

NSUME 

SUMA 


Asimuth (degrees) input data. Output value placed 
in AZ(5) . 

Elevation (degrees) input data. Output value 
placed in EL(5). 

AZ rounded to nearest degree. 

EL rounded to nearest degree. 

DO- loop index. 

Index value of the mode (highest population) inter- 
val, azimuth. 

Index value of the mode (highest population) inter- 
val, elevation. 

DO-loop index. 

DO- loop index. 

"Distribution" of lAZ, number of AZ values rounding 
to the correspondingly indexed lAZ (includes harm- 
less extraneous values) . 

"Distribution" of lEL, number of EL values rounding 
to the correspondingly indexed lEL (includes harm- 
less extraneous values) . 

Population of the five-degree azimuth averaging 
interval centered on the mode lAZ(KMA). 

Population of the five-degree elevation averaging 
interval centered on the mode lAZ(KME). 

Sum of the AZ lying in the averaging interval for 
azimuth. 
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SUME 


Sum of the EL lying in the averaging interval for 
elevation. 
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SUBROUTINE TRACK 


De>cription 

The purpose of Condenser (MAIN) is to discern from the O.I- 
second raw data, first, the switch points, i.e., the points at 
which the baroswitch changes contacts (channels), and second, a 
condensed representation of the signal transmitted while on each 
contact. The representation consists of the three quantities for a 
given signal dwell: 

CONDd, JK) * switch, or beginning, time (seconds from launch) 
C0ND(2, JR) > duration (seconds) 

CONDO, JK) ■ mean frequency 

If the signal changes suddenly for a given contact position, more 
than one signal dwell may result for that contact position. Thus, 
a rapidly varying or noisy signal may be represented by several 
condensed points, COND, over its baroswitch interval. 

Subroutine TRACK examines the ten-point sample, FREO(i), i ■ 
1, 10, to determine whether the signal is in the frequency tracking 
gate; i.e., whether at least two of the FREQ(i) lie within HGATE 
(Hz) of SIGLEV (Hz). If so, SIGLEV is adjusted half-way toward the 
mean value of the points within the gate. This mean value includes 
SIGLEV for further stabilization. Some computing time is saved by 
requiring no more than the first seven points in the gate for com- 
putation of the mean. TRACK continues by returning to ADVANC via 
MAIN (without changing JK, i.e., JK • JKMEM) for more data and 
repeating the process, summing the mean values: 
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I 



FSUM ■ FSUM * SUM6TE/NGATE 
MSUM -NSUM -t- 1 

uatil the signal is not found in the gate. 

If less than two of the ten-point raw data samples lie in the 
frequency tracking gate, TRACK calls subroutine SEARCH. SEARCH 
obtains successive ten-point samples, FREQ(i), by returning to 
ADVANC via TRACK and MAIN, without changing JK (JK > JKMEN) . If 
SEARCH cannot find signal or exceeds COND dimension, JK > 1000, it 
returns to TRACK (S.83, S.85, respectively), and TRACK returns to 
Terminator (MAIN S.83, S.85, respectively). 

When SEARCH finds signal, it returns to MAIN via TRACK, hav- 
ing incremented JK (JK * JKMEM) , where control therefore proceeds 
to subroutine DCCOM to process the new condensed point, C0ND( , 
JK). 







I 
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Plow Diagram, TRACK 






CALL Lilt, TRACK (Ref.! List of Voriablet, below) 


Variable 

Flow 

Coomenta 

TIHE(IO) 

From ADVANC via MAIN 

Ten-point raw data sample, 
TIME sent to SEARCH, FREQ, 

FREQ( 10) 

From ADVANC via MAIN 

used to track signal. 

TEST(IO) 

From MAIN (Initializer) 

Diagnostic printout 
control, input TEST(i), i ■ 
1-3 used. 

LOS 

To MAIN (Terminator) 
from SEARCH 

Counts continuous half- 
seconds of no signal. 

CONDO, 1000) 

To DECOM via MAIN 
from SEARCH 

Condensed points generated. 

JK 

Variables in 
COMMON/ SIGNAL/ 

To MAIN 

« 

• 

Index of last computed 
COND. 

SIGMAX, 
SIGMIN, 
HGATE, IN 

From MAIN (Initializer) 

For SEARCH 

SIGLEV 

From SEARCH, initially 
from MAIN (Initializer) 

Frequency gate position 

NSUM, FSUM 

From/To SEARCH 

For mean frequency computa- 
tion. 


- 44 - 


Li«t of Variablet, TRAOC 

BLWR, BUPR Lower end upper bounds (Hz), respectively, of the 
signel trecking gete. 

C0ND(3, 1000) Array (real) of condensed points defined in SEARCH. 

FREQ(IO) Ten-point sample of the O.l-second raw data (Hz), 

used to determine whether the signal is in the 
gate, and, if so, to center the gate and to compute 
SUMGTE and FSUN. 

FSUM Running am over a dwell of the mean frequency in 

the signal tracking gate, used by SEARCH in comput- 
ing CONDO, ). 

HGATE Half-width (Hz) of the signal tracking gate. 

IN The number of subdivisions of the sonde signal 

frequency range, used by SEARCH. 

J DO-loop index. 

JK Count, index of the last condensed point computed, 

incremented by SEARCH, used by MAIN, etc. 

LOS Counter of consecutive half-second steps for which 

no signal is detected by SEARCH. 

NGATE The number of the ten raw data FREQ which fall 

within the signal tracking gate, SIGLEV t HGATE, 
augmented by one to include SIGLEV in the computed 
mean. 

NSUN The nisnber of half-second samples in a signal dwell 

which %fere found in the gate, used in computing the 
mean frequency CONDO, ) ■ FSUN/NSUM. 
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SIGLBV 


SIGMAX, 

SIGMIN 


SUMGTE 


TEST(IO) 

TIME(IO) 


Center, poeition (He), of the eignel Cracking gate, 
eucceesively adjusted with lag to the signal mean; 
i.e., SIGLEV - (SIGLBV * SUHGTE/NSUM)/2. 

Upper and lower limit (He), respectively, of the 
sonde signal frequency range, used in excluding 
noise. (SIC^IAX > 205, SIGMIN -5) 

Sun of the frequencies FREQ which lie within the 
signal tracking gate, plus the preceding value of 
SIGLEV, used to compute local means SUMGTE/NGATE 
which are again averaged over the dwell, %ihich are 
used to adjust the gate position SIGLEV at each 
step. 

Diagnostic printout control, input, TRACK uses 
TEST(i), i » 1-3. 

Times (seconds from launch) corresponding to the 
ten-point raw data sample FREQ, used by SEARCH. 


46 - 


SUBROUTINE SEARCH 


De»cription 

When a sudden change in signal level occurs, i.e., when a 
ten-point raw data sample falls outside of the signal tracking gate 
SIGLEV ± HGATE, subroutine TRACK calls subroutine SEARCH to find 
the new signal level, and to *'wrap up" the condensed data point 
just passed. SEARCH searches the signal range from SIGMIN to 
SIGMAX to find the frequency band of half-width HGATE which con- 
tains most but at least three of the ten points (Fig. 3) . If more 
than one of the bands contain the largest number of points, the one 
at highest frequency is chosen. 

If no band contains at least three points, LOS is increment- 
ed, and SEARCH continues by returning to ADVANC via TRACK and MAIN 
for more data. This proc^^ss is repeated, except: 

a. At LOS > 1, TSWCH2 (the beginning time of the present 
noise or LOS interval) is assigned (this quantity is used 
as the ending time of the preceding signal dwell, thus 
excluding the noisy interval from that signal dwell), and 

b. At LOS > 100 (30 seconds), indicating excessive noise or 
no signal, at %diich time SEARCH terminates processing by 
returning to TRACK and MAIN (S.83). 

If, as in the normal case, a frequency band is found suffi- 
ciently populated, the center (SIGLEV) of the signal frequency gate 
is placed at the midpoint of the band, the beginning time of the 
new signal dwell is assigned, LOS is stored in LOSN (for use in 




- 47 


Index IB 


Frequency 


IM - — ^ - SXGMAX ■ 2.5 Hz (*100 ordlnetea) 






3 

2 

1 



HGATE ■ 1.0 Hz (*0.5 ordinates) 


SIGMIH ■ 5 Hz (*2.5 ordinates) 


Index 


Count 


Frequency (Hz) 


■■ ( 

KBND<^ 

1 

KB 


KBL 



KBLL 






(b) 


Fig. 3. "Met" signal acquisition, (a) To find signal, the entire 
sonde frequency range is divided into IN one*hertz 
levels, (b) Beginning at SICMIN ♦ 2 * HGATE, the nvssber 
KB of the ten raw data points FREQ falling below each 
level BHD, less that KBLL of two levels below, is cois- 
puted. KBND ■ KB - KBLL. The index IBND and count KBNDC 
of the most populated such 2-Hz band is obtained. If 
KBNDC is at least three points, the signal-tracking gate 
SIGLEV is centered at the level corresponding to index 
IBND - 1; i.e., 

SIGLEV - SIGNIN * (IBND • 2) * HGATE 
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computing DWELL) before being reeet Co tero, COND index JX is in* 
cremented, and a new condensed point is defined over the preceding 
signal dwell. 


CONDd, JK) • TSWCHl 
C0N0(2, JR) - DWELL 
CONDO, JR) - FSim/NSUM 


Dwell 

Noise 


or L.O.S. 


TSWCHl ;SWCH2 TSWTCH 

TIME(l) 


The signal frequency over the dwell is taken as the mean of the 
’’ten-point” means FSUM computed in TRACK. 
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Plow Di«gr— . SEARCH 


' JLW-. 


1 



Ooternine 2 Hz frequency band of signal, IBHD, 
and nuaber of points, KBNOG, in that band 



1 


i 



SO 















CALL List, SEARCH (Ref.; List of Variables, below) 


Variable 

Flow 

Comments 


TIME(IO) 

From ADVANC via 
MAIN, TRACK 

Ten-point 

raw data sample 

FREQ( 10) 

From ADVANC via 
MAIN, TRACK 

Ten-point 

raw data sample 

LOS 

To MAIN via TRACK 

Count of 
seconds of 

continuous half' 
no signal 

C0ND(3, 1000) 

To DECOM via TRACK, 
MAIN 

Condensed 

points generated 

JK 

To MAIN (Terminator) 
via TRACK 

Index of last COND computed 

TEST(IO) 

From MAIN (initializer) 
via TRACK 

Diagnostic 

input 

printout control 


In COMMON/ SIGNAL/ 


SIGMAX, SIGMIN, 

HGATE, IN From MAIN (Initializer) Signal range and increments 
via TRACK for search 


SIGLEV 


To TRACK (initially For TRACK 
from MAIN via TRACK) 


NSUM, FSUM To/From TRACK 


Accumulators for mean fre- 
quency computation 
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List of Variables , SEARCH 


BND Momentary upper bound (Hz) of a signal search re- 

gion. 

CONDO, 1000) Array (real) of condensed points defined by Con- 
denser (TRACK and SEARCH), each being the dwell or 
interval between sudden changes in the signal: 
CONOd, ) ■ beginning time (seconds from launch) of the 
signal dwell. 

* duration (seconds) of the signal dwell. 

* mean signal frequency (Hz) over the dwell. 
Duration (seconds) of the current signal dwell. 
Ten-point raw data sample, used to locate the sig- 
nal in the search process. 

Sum over the dwell of the mean frequency values 
within the tracking gate computed by TRACK, used to 
compute the mean frequency C0ND(3, ). 

Step size (Hz) used in the search for signal over 
the sonde signal range. 

DO-loop index, used in the stepping search process. 
Index of the most populated frequency band when 
searching for signal, each bandwidth is two steps 
HGATE. 


C0ND(2, ) 

C0ND(3, ) 

DWELL 

FREQ(IO) 

FSUM 


HGATE 

IB 

I BND 


IN The number of subdivisions (steps) of the sonde 

frequency range used in searching for signal, IN ■■ 
((SIGMAX - SIGMIN) /HGATE) + 1, defined in MAIN 
( Initializer) . 


- 52 - 


ITYPE 

J 

JK 

KB, KBL 

KBND 

KBNDG 

LOS 

LOSU 

NSUM 

SIGLEV 

SIGMAX, 

TEST(IO) 

TIME(IO) 

TSWCHl 


Program mode (see Input Card Deck). 

DO-loop index. 

COND index, array counter. 

KBLL Counter of raw data FREQ below level BND in the 
signal search process (see SEARCH, Description). 

The number, 1-10, of raw data FREQ falling within 
the 2*HGATE interval below BND (Hz). 

Largest KBND over the 5-205 Hz range of the sonde. 
Counter of consecutive half-second steps for which 
no signal is detected by SEARCH. 

Stored value of LOS, used to recall LOS > 0 to 
exclude noise interval from DWELL. 

Accumulated number of half-second raw data samples, 
in the signal dwell, found in the tracking gate by 
TRACK. 

Tracking gate center (Hz), repositioned on the 
signal by SEARCH. 

SIGMIN Upper, lower limits (Hz), respectively, of the 
sonde signal frequency range. 

Diagnostic printout control, input, SEARCH uses 
TEST(i), i -1-5. 

Time (seconds from launch) corresponding to the 
ten-point raw data samole FREQ, used to identify 
sudden signal changes; i.e., dwell boundaries. 
Beginning time (seconds from launch) of the current 
dwell. 
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TSWCH2 


TSWTCH 


Beginning time (seconds from launch) of a noise (no 
signal) interval, used as the ending time of the 
dwell when succeeded by noise. 

Ending time (seconds from launch) of the current 
dwell. 
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SUBROUTINE DECOM 


Description 


Subroutine DECOM determines for each condensed point, COND, 
its proper channel (Temperature, Reference, High Reference or 
Humidity) and for each Reference and Humidity switch point, its 
baroswitch contact msaber. Channel and contact numbers are stored 
in ICOND. 


lCOND(i, ) ■ 1 Temperature 
■ 2 Reference 
* 3 High Reference 
» 4 Humidity 

> 5 Undecommutated, rejected 
1C0ND(2, ) “ 1-180 Baroswitch contact number 

Decom may be considered in two major parts, Process "A" in 
which Temperature and Humidity channel and contact numbers are 
assigned as each occurs ("forward-assigned”) (S.178, S.200), and 
Process "B" in which, after assigning each Reference its channel 
and contact number (S.600), the preceding Temperature and Humidity 
channel and contact numbers are corrected ("back-assigned") where 
necessary (S.180). Process ”B” also discriminates High Reference 
dwells and verifies their contact numbers (S.130), and detects 
balloon burst (S.30). These parts are indicated in Fig. 4. More 
detailed discussion follows. 


- 55 - 


1 ^ 



Process "A 


Process "B 


Overall block diagram of subroutine DECOM. 







Initialise 


Initial values for TNOH, DSL. FHUM, ICRl in the CALL Liat are 
set in MAIN. TNOH ■ 10000.0 is set large so that it is ineffective 
as a burst indicator in INTER? unless TNOH is computed in OECOM. 
Starting frequencies DSL and FHUM for the Temperature and Humidity 
tracking gates are computed in MAIN according to input ordinate 
values at launch. The effective intercontact value ICRl of the 
baroswitch at launch is computed in MAIN according to the input 
surface pressure FP0. 

Initial values of RFL and DRPFL, e.g. , 170. and 2. Hz, have 
proved acceptable. PFL is set initially somewhat large to prevent 
its premature interference, and low enough that it reaches a proper 
value by the time it is needed at contact 30. Accumulators RFSUM 
and RTSUM and reference marker TR are initialized zero. Contact 
number multiplier MLT and its real form AMLT are set at the con- 
stant 100. NXTP is initially made 30, the contact number of the 
first high reference. 

Temperature tracking gate half-width, GTEMP, is initially 4.0 
Hz (GTEMP is subsequently increased to 6.0 above contact 135 where 
Humidity data ceases, allowing greater changes in lapse rate be- 
tween the larger Reference intervals later in flight). Ml is made 
unity as if under normal conditions. Quantities TSL, ESL, TRl, 
INCH, JKR, NOH, JKP, ICM, JRRl, ICR, and KROSS are initialized 
zero. 

Upon each call of DECOM, the variables T, D, DWELL, D, and TB 
are defined from the condensed point COND to be processed. 
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Forward-Aiiign 


Channel Pecommutation . 

Contiguous condenaed points or dwells whose frequencies ex- 
ceed the reference threshold RFL are accumulated (S.220) for a mean 
frequency and are assigned channel nisnber 2. The channel number 
for each condensed point is also remembered in INCH, for use in 
processing the next condensed point. The index of the first dwell 
of a Reference group is retained in JRR (S.200). Upon arrival of a 
reference DWELL, C0ND(3, ) ■ D > RFL, the value INCH * 2 indicates 

that the arriving dwell is a continuation of the preceding dwell, 
whereas INCH 0 signifies that the arriving dwell is the first 
dwell after launch. If INCH is neither 0 nor 2, the arriving dwell 
is taken as the beginning of a new Reference group of dwells; i.e., 
a baroswitch contact switch point. The group may consist of more 
than one dwell (COND element), depending on the variation of the 
signal over the group. 

Signal dwells which are not Reference points, C0ND(3, > D < 

RFL, which immediately follow a Reference, INCH ■ 2, if not after 
the first Reference of the flight, JK0 ¥ 0, send control to Process 
"B" (S.600) discussed below. If after the first Reference, JK0 • 0, 
the point is taken as a Temperature dwell. Points assigned to the 
Temperature channel are used to adjust the temperature tracker 
variables ESL (S.7S2), DSL (S.781), and TSL. The temperature 
tracking gate position at a subsequent time TB is computed (S.78): 

TF - DSL ♦ F.SL*(TB - TSL) 
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Dwell frequencies D felling within GTEMP of TF ere eccepted es 
Tefflpereture dete. Dwells felling outside this gete but inmedietely 
following e Reference, i.e., when 


NOH - 0, JK0 00, T < TR + 0.8*SLOPE 


or when 


NOH - 1, INCH » 2 

ere accepted as Temperatures and are used to correct or "recapture" 
the temperature gate (S.7831). NOH “ 0, 1 signifies the humidity- 
no humidity regions of flight defined by contact number 135. Dwells 
not accepted under the above conditions are processed as Humidity 
data (S.785) if NOH ■ 0, or are rejected (S.784) if NOH « 1. 

Inter-Reference Contact Number Determination 

Dwells to be processed as Humidity (S.785) which are not 
"switch points", INCH * 1, are accepted as "data points" of Humidi- 
ty groups (S.788). If, however, the dwell follows a Temperature, 
INCH * I, it is treated also as a baroswitch contact switch 
point. If it is the first such point, ICM ■ 0, the contact number 
is computed as the next integer greater than the decimal contact 
number, ICRl (ICR0 in MAIN), associated with the surface pressure 
at balloon launch. The contact time rate SLOP2 is also computed 
for use in computing subsequent contact numbers. Subsequent con- 
tact numbers are computed according to elapsed time, (T-T2), from 
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the preceding contact ninber ICM # 0 and according to the contact 
time rate 8LOP2. After t%m Reference eifitch pointa have occurred, 
SL0P2 ie computed from the preceding inter-Re ference time inter- 
val. When a computed Humidity contact mnber exceeds the next 
expected Reference contact number, the switch point is rejected, 
ICONDCl, ) ■ 6, and "back-assigning" is requested, KROSS ■ 1, over 
the entire inter-Re ference interval. Such a condition is caused by 
erroneous decommutation (channel crossovers), by momentary balloon 
descents ("dip") , or by baroswitch or signal irregularities. 

Process Reference 

When a non-Re ference dwell is encountered immediately fol lott- 
ing a complete (JiU) ^ 0) Reference group (INCH ■ 2), then control 
is sent to Process "B" (S.600) where, first, the Reference group is 
processed. The mean frequency and cumulative dwell time over the 
Reference group are assigned to the leading dwell C0ND( , JKR) . 

If it is the first Reference group (ICR * 0), its contact 
number is computed as the next multiple of five integer greater 
than the decimal contact number at launch (ICRI), plus an addi- 
tional multiple of five for each whole multiple of 90 seconds 
contained in the time interval from launch to TR seconds. 

ICR » (ICRBl ♦ ICRB) * 5 * ML 


where 
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ICRBl ■ ICR1/(5*HLT) 

ICRB ■ TR/90. +0.5 
> 1 

Ninety seconda is a sufficiently valid estimate for the inter**Ref* 
erence time interval. The resulting contact number ICR accommo- 
dates missing first References. 

Succeeding Reference contact nimbers are incremented by 5*M 
for contact msnbers less than or equal to 135 and by H for contact 
numbers greater than 135. The intercontact time interval SLOPE is 
computed assuming no Reference switch points were missed (M ■ 1) 
and compared to the previous value (SLOPl) to verify or determine 
the correct value for M. 

M - SLOPE/ SLOPl +0.5 

When M P 1 and recomputation of ICR is required, the fact is remem- 
bered, Ml N P 1, to inhibit later "back-assigning" over such 
inter-Reference intervals. M-1 is the nvmber of References miss- 
ing. 

When M < 1, signifying too early a Reference contact, baro- 
switch reversal is assumed (by setting Ml ■ -1) and the switch 
point is ignored. Such cases can be due to balloon "dip" 
(temporary descent). 

After accepting a Reference switch point (S.645), GTSW is 
computed for use in back-assigning, described below. At contact 
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number 135 (ICR ■ 13500) the NOH flag it set to unity, the 
corretponding tine TNOH set, and the temperature gate half-width 
GTEMP it increated to 6 Ha. 

If the mean frequency of the Reference group, C0ND(3, JKR) 
(S.120), it lett than the High Reference threahold, PPL, both PPL 
and RPL are adjuated to follow: 

RPL - 0.6*RPL ♦ 0.4*(COND(3, JKR) - 10.0) 

PPL - 0.6*PPL ♦ 0.4*(RFL ♦ DRPFL ♦ 10.0) 

Otherwiae the Reference it taken at a High Reference (S.130) and 
the Reference-High Reference difference it adjuated: 

DRPPL - 0.6*DRPPL + 0.4*(C0NO(3, JKR) - RPL - 10.0) /2 

If the High Reference is not the first one, tests are made to de- 
termine whether a High Reference was missed, and the next expected 
High Reference contact mnber NXTP it defined. 

Finally, in processing Reference points, the dwell, CONlK , 
JKR-1), immediately preceding the Reference is again defined a 
temperature datum. If this requires a change in ICONOCl, JKR-1), 
then the temperature frequency gate it "recaptured" and back- 
assigning it requested (RROSS * 1). 

Test for Burst 


When Reference switch time TR exceeds 3000 seconds (50 min- 


utes from launch), *,\to tests for bslloon burst are performed. The 
principal test (S.9S1 * 1) simply computes the time (TBRST) corre- 
sponding to the manually input decimal contact number at burst 
(CBRST). That is, when ICR reaches the contact number just preced- 
ing burst, TBRST is computed according to the fraction of the con- 
tact interval remaining before burst. 

TBRST ■ TR ♦ (CBRST - AINT( CBRST) )*SLOPE 

The alternate test automatically determines the '.ime of burst 
by detecting the sudden decrease in intercontact time SLOPE. A 
running mean (SDWELT) over the preceding four contact periods 
(SLOPE), excluding SLOPE < 15 seconds (interpreted as barocwitch 
noise), is examined at each reference switch time to detect a 
sudden decrease. When three successive decreases (IBRST “ 3) in 
this smoothed SLOPE occur, burst time is assigned to the third 
preceding contact switch point. TBRST « T0WELT(4) 


|— — SDWELT (4) (3) 

TDWELT (4) 

t 

TBRST 


“4— (2) —4^ (l)-*j 

(3) (2) (1) 


SLOPE decreasing, 

■ ■ — ' - ' ' 

sonde falling 


Manual burst input CBRST routinely should be supplied to terminate 
data processing before burst to avoid spurious data conditions at 
apogee. 
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B«ck-Atiign 

When indicationt of error occur in the forward-aasiining 
process above, back-assigning is automatically requeated (KROSS ■ 
1) for the respective inter-Re ference interval. Back-assigning 
incorporates the additional advantages of using the value of 
contact time rate, SLOPE, computed over the inter-Re ference 
interval to %#hich it is being applied (instead of using SLOPE 
computed from the preceding interval) , and of using time gates as 
%fell as frequency gates for diacriminating Temperature and Humidity 
points. Back-assigning is ineffective, however, if the contact 

rate is irregular as in the cases of balloon dip or missed 
reference switch points. In such cases back-assigning is prevented 
(Ml * 1). 

The nearest contact number, NCT, to a given dwell is com- 
puted, along with its estimated time, Tt. 

TS • (TN - TRl)/ SLOPE 
NCT • TS ♦ 0.5 
TO • FLOAT(NCT)*SLOPE ♦ TRl 

The expected frequency, TFN, of the temperature signal is computed 
from the temperature tracking parameters DSLN and CSLN computed 
previously st TSLN, 

TFN • DSLN ♦ ESl.N»(TBN - TSLN) 



! 






A dwell which falls in the frequency gate, 


TFN ± GTEMP 


or lies on either side of a Reference (MCT " 0, 1, 5), and whose 
midtime TEN falls GTSW seconds before its nearest expected contact 
switch time T0, is accepted as a temperature point. Before 
proceeding to the next older dwell, the temperature tracking param- 
eters are updated according to the frequency DN and time TEN of the 
accepted temperature point. 



I 


ESLN » 0.8*ESLN + 0.2*(DN - DSLN)/(TEN - TSLN) 

DSLN » 0.8*DSLN + 0.2*DN 
TSLN * TEN 

except if ESLN changes too suddenly, more than 0.2 Hz/s, it is left 
unchanged. This protects against gate-stealing by noise. 

A dwell not accepted as a temperature point is next consid- 
ered as a contact switch point. If: 

a. Its leading edge tine TN falls within GTSW seconds of its 
nearest expected contact switch time T0, and 

b. Its frequency is not within the temperature frequency 
gate, and 

c. Its nearest expected contact switch time is between (not 
including) the References (0 < NCT < 5), and 

d. It is the l.eading dwell satisfying these conditions for a 
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given contact nimber NCT, 

then it is accepted as a contact switch point. It is assigned 
contact number 


1C0ND(2, ) « 1C0ND(2, JKRl) ♦ NCT 

and channel number 4 (Humidity) datum. 

A dwell not accepted as a temperature nor a switch point 
under the conditions above, but was accepted as a temperature in 
the forward-assigning process, is accepted as a temperature 
point. All other dwells treated in back-assigning are rejected. 

Before resuming the forward-assigning process (S.78), the 
quantities JKRl, TRl, and JNSTRT are advanced in case back-assign- 
ing is requested for the next inter-Re ference interval, ICRl, 
SLOPl, T2, SL0P2 are advanced for use in forward-assigning, and Ml, 
KROSS are reinitialized (S.188). 

Optional Diagnostic Print 

Printout during execution of subroutine DECOM includes op- 
tional as well as warning messages. When input TEST(9) is unity 
(greater than 0.01), internal back-assign quantities are printed. 
Until dwell time T exceeds input TEST(8), forward-assign quantities 
are printed at each exit of DECOM. 
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(Process "A") (continued) 

























DECOM (Process "B”) (continued) 
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IC0ND(2,JN) = ICOND(2,JRKl) + NCT 
















CALL List, DECOM (Ref.: List of Variables, below) 


Variable 

CONDO, 

JK 

IC0ND(2, 

TNOH 

DSL 

TEST(IO) 

ICRl 


Flow 


1000) From TRACK via MAIN 

From SEARCH via MAIN 
and TRACK 

1000) To INTER? via MAIN 


To MAIN for INTERP 
(large initial value 
from MAIN) 

Initial value from MAIN 


From MAIN (card input) 


From MAIN (initial 
value) 


Comments 

Condensed data points 
Index of C0ND( , JK) 


Channel and contact num- 
bers of COND, computed by 
DECOM 

Time at baroswitch 
contact number 135 


Temperature gate position 
at launch 

Controls for diagnostic 
printout 

Effective contact nvmber 
at launch 
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List of VTi«ble>, DECOM 


CBRST 

COND(3, 1000) 
CONDd, ) 


AMLT Multiplier, reel form (AMLT “ 100.), used in the 

computation of baroswitch contact numbers for more 
precision. 

Effective baroswitch contact manber at balloon 
burst, real form (XXX. XX). 

Real array, condensed data: 

■ elapsed time from launch (seconds) to the leading 
edge of the dwell. 

C0ND(2, ) ■ duration (seconds) of the dwell. 

CONDO, ) •• mean signal frequency (hertz) of the dwell. 

D, DN Same as C0ND(3, ), used in forward- or back-as- 

signing processes, respectively. 

DRPFL Running difference (Hz) bet%raen High Reference and 

previous low Reference. Initial value set equal 2. 

DSL, DSLN Temperature gate position (Hz) computed and used in 

forward- or back-assigning processes, respectively. 

DWELL, DWELLN Same as C0ND(2, ), used in forward- or back-as- 

signing processes, respectively. 

DWELT Baroswitch period (seconds/contact), used in auto- 

matic detection of balloon burst. 

ESL, ESLN Slope (Hz/second) of Temperature signal frequency, 

used in first-order extrapolation of Temperature 
gate, in forward- or back-assigning processer, 
respectively. 

ESLl, ESLNl Previous value of ESL, used if computed value ex- 
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ceed* a maximum change, in forward- or back-aatign- 
ing procesa. 

GTEMP Half-width (Ha) of Che Temperature gate. 

GTSH Half-width (aeconda) of the contact awitch gate, 

uaed in back-aasigning procesa. 

1 Index (arbitrary), uaed in printing statement. 

IBRST Counter used in automatic detection of burst. 

ICM Integer indicating baroawitch contact manber, mul- 

tiplied by MLT ■ 100, 0 < ICM < 18000. 

IC0ND(2, 1000) Integer array, condensed data: 

ICONDd, ) ■■ channel mmber: 

1 Temperature 

2 Reference 

3 High Reference 

4 Humidity 

> 5 Undecommutated 

IC0N0(2, ) ■ baroawitch contact number, 1-180. 

ICR Computed contact number of the current Reference 

point, multiplied by MLT > 100. 

ICRl Stored value of preceding ICR. 

ICRB Number of Reference contacts from ICRBl to the 

first one detected. 

ICRBl Contact number of the Reference contact "lower" 

than that at launch. 

IDC Change in contact mraber since the last detected 

contact switch time. Usually unity. 
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INCH 
IS, ISS 

ITCNT 

JK 

JKP 

JKR 

JKRl 

JN 

JNl 

JNFIN 

JNLCT 

JNSTRT 

KROSS 


LCT 


M 


Ml 


Channel of preceding dwell. 

Indices used in automatic detection of balloon 
burst. 

Counter used for labeling diagnostic printout, 
TEST(8). 

Index of condensed point, or dwell, COND( , JK) , 
being processed by OCCOM (in forward-assigning). 
Index of the previous High Reference switch point. 
Index of the current Reference switch point. 

Index of the previous Reference switch point. 

Index of condensed point, or dwell, COND( , JN) , 
being processed in back-assigning process. 
JN-decrementing variable in back-assigning DO-loop. 
First, largest, index JN in back-assigning process. 
Value of index JN at preceding humidity point. 

Last, smallest, index JN in back-assigning process. 
Flag requesting (KROSS ■ 1) back-assigning. Set 
under conditions which indicate errors in forward- 
assigning process, e.g., T-H crossover, etc. 

Value of contact number IC0ND(2, JNLCT) assigned to 
the preceding humidity point. 

Integer used in computing contact number, repre- 
sents increment of contacts according to elapsed 
time. 

Flag indicating irregular contact progression 
(Ml * 1) which contraindicates bad ^assigning. 
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MLT Multiplier constant (MLT ■ 100) which serves to 

increase precision in contact nunber computation. 

MM Inert tag used for diagnostic purposes in back** 

assigning process. 

NCT Integer used in back^assigning process. Represents 

contact number increment from the earlier Reference 
contact number lC0ND(2t JKRl). 

NOH Integer indicating state below (NOH * 0) or above 

(NOH 1) contact number 135. 

NXTP Number of next expected High Reference. 

PFL High Reference threshold (Hz), used to discriminate 

High and ]ow Reference. 

RFL Reference threshold (Hz), used to discriminate 

Reference dwells from Temperature and Humidity 
dwells. 

RFSUM Summation variable (Hz) for computing the mean 

frequency over a Reference group. 

RTSUN Summation variable (seconds) for computing the mean 

frequency over a Reference group. 

SOWELT Running mean over four DWELT. 

SLOPl Stored SLOPE from preceding inter-Reference inter- 

val, used in computing (forward-assigning) Refer- 
ence switch point contact number. 

SL0P2 Stored SLOPE from preceding contact interval, used 

in computing ( forwanl-assigning) Humidity switch 
point contact number. 
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SLOPE 


T, TN 


T0 


T2 


TB, TBN 


IBM 


TBRST 


TOWELT 


TEST( 10) 


TF, TFN 


Mean contact time rate (aeconda/contact) between 
the preceding two Reference awitch points. 

Switch time, time (seconds from launch) of leading 
edge of the current dwell for condensed point. 
Same as C0ND(1, ), in forward- or back-assigning 

procei*ses, respectively. 

Switch time of nearest expected contact, in back- 
assigning. 

Switch time of preceding contact, in forward-as- 
signing, used to compute contact number from 
elapsed time. 

Midpoint time C'T-bar”) of the current dwell (sec- 
onds from launch) . 

Time of burst (minutes from launch) from automatic 
detection of balloon burst, for auxiliary printout 
only. 

Time of balloon burst (seconds from launch) com- 
puted from CBRST or from automatic detection of 
burst. 

Time (seconds from launch) corresponding to WELT 
and SDWELT. 

Input control of diagnostic printout, TRST(6) and 
TEST(9) only, used in DECOM. 

Expected frequency (Hx) of the next Temperature 
dwell, center frequency of Temperature gate in 
forward- cr back-assigning processes, respectively. 


TNOH 


1 

i 

TiiM (••condt fron launch) of contact 115, at which 
humidity data terminataa. 

TPl Time (aeconda from launch) of the preceding High 

Reference awitch point, used in detecting and cor- 
recting for missed High Reference points. 

TP12 Expected time (seconds) between High References, 

used in detecting and correcting for missed High 
References. 

TR Switch time (seconds from launch) of the current 

Reference point. 

TRl Switch time (seconds from launch) of the preceding 

Reference point. 

TS Estimated number (reiil) of contact numbers from TRl 

to TN, used in computing the nearest contact nunbex 
in back-assigning. 

T5L, TSLN Time (seconds from launch) of the previously-com- 

puted Temperature dwell from which the Temperature 
gate is extrapolated, used in forward- or back- 
assigning, respectively. 
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SUBROUTINE INTER? 


Description 

The primary purpose of INTER? is to construct a table of 
Pressure (mb). Reference frequency (Hz), Temperature (ordinates), 
and Relative Humidity (ordinates) values at one-minute intervals 
from the surface to balloon apogee (burst), to be used as input 
data to ECC-PRD. NASA computer program ECC-PRD, used without 
ozonesonde inputs, produces the final output of the RAWINSONDE data 
processing system. 

INTER? receives decommutated sonde data as asynchronous 
samples COND( , i), ICOND( , i), i = 1 to JK, the one-minute table 
to be completed VL(i, j) , i “ 4 to 7, j * 1 to LIST, surface values 
V2(i), i “ 4 to 7, the baroswitch pressure calibration function 
PCAL(i), i * 1 to LCNTK, the time of flight at the end of humidity 
data, TNOH, and at the end of all data, TBRST, and the diagnostic 
printout control input, 'TEST(6). INTER? computes the values which 
complete the table VL and defines ISTOP “ 10 when it reaches the 
time of balloon burst (TBRST) in the data, or ISTO? ® 2 if it 
reaches the last baroswitch contact calibrated (LCNTK). The last 
contact number used by INTER? is sent back to MAIN for printout in 
case ISTOP “ 2. 

Interpolation 

Linear interpolation is performed in each of the four vari- 
ables: contact number, Reference (Hz), Temperature (ordinates), 

and Humidity (ordinates). Bracketing values TI(IV), VI(IV), 

- 79 - 


i 

j 


I 


i 


I 


i 

i 



T2(IV), V2(lV),in the condensed data (COND, ICOND) for a given one- 
minute level, L, are accepted for each variable according to the 
following tests. 

For contact mmber (IV ■ 4) only the first-occurring time of 
a given contact number is used, and no contact number is used if 
the condensed point was finally determined undecommutated, 
ICONDd, ) >5, or if the contact number exceeds the highest con- 
tact mmber calibrated, LCNTK. 

For Reference frequency (IV » 5) , only the frequency of Ref- 
erence switch points (those whose frequencies were computed by 
DECOM over the entire group of reference dwells) are accepted. 

For Temperature (IV ■ 6), only dwells which have nonzero mean 
frequency are accepted. The frequency value is converted to ordi- 
nates with the use of the local Reference frequency VL(5, L) ac- 
cording to 


V2(6) “ 95.*COND(3, l)/VL(5, L) 

The time is taken as that at the midpoint of the dwell. 

T2(6) = CONDd, I) + C0ND(2, I)/2 

For levels beyond the last temperature dwell, e.g., during a Refer- 
ence dwell just before balloon burst, the extrapolating quantities 
are chosen to be the one-minute values at the two preceding 
levels. This is done to avoid possible large errors resulting from 
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extrapolation from ahort-* tine-base variable temperatures which may 
occur in the data. 

For Humidity, frequencies less than 5 Hz are excluded. The 
time T2(7) and ordinate V2(7) are computed in the same way as for 
Temperature. No values are computed beyond contact 135, i.e., 
beyond time TNOH. 

For all four variables, no condensed data (COND, ICONO) is 
used which exceeds JK in index or TBRST in time. 

Diagnostic Printout 

Interpolation quantities during the first and last 30 minutes 
of the flight are printed for diagnostic purposes if input value 
TEST(6) is greater than 0.01 (e.g., TEST(6) » 1). Included in this 
printout are the input quantities JK, LIST, TNOH, V2, and TBRST. 

Signal Dropout 

Zeroes are entered in the one-minute table (VL) at levels 
’.here no measurement data, C0ND( , ), are within ALOSS seconds of 

the one-minute level. ALOSS is 200., 600., 100., and 100. seconds 
for pressure. Reference frequency. Temperature, and Humidity, re- 
spectively. 














INTER? (continued) 























CALL Li>t, INTER? (Ref., List of Variables, below) 


Variable 

Flow 

Comments 

1C0ND(2, 1000) 

From DECOM via MAIN 

Channel, contact number, 
asynchronous data 

C0ND(3, 1000) 

From Condenser, SEARCH 
via TRACK and MAIN 

Time, dwell, frequency, 
asynchronous data 

JK 

From condenser, SEARCH 
via TRACK and MAIN 

Count, of COND, ICOND, 
condensed points 

PCAL( 180) 

From Initializer MAIN 

Baroswitch contact cali- 
brated pressures 

TNOH 

From DECOM via MAIN 

Time at contact 135, 
when humidity data stops 

TBRST 

From DECOM via MAIN 

Time of burst, apogee 

ISTOP 

To Terminator MAIN 

Terminal condition in- 
dicator 

LCNTK 

From Initializer MAIN 

Highest contact number 
calibrated 

KNTCT 

To Terminator MAIN 

Last contact number used 
by INTERP 

V2(7) 

From Initializer MAIN 

Upper bracketing quanti- 
ties, become initial 
(surface) values of VL 

TEST(IO) 

From Initializer MAIN 

Diagnostic printout 
control, input 
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List of 
AL0SS(7) 

CONDO, 

DUST 

1 

14. 15, 

1C0ND(2, 

IJ 

lOUT 
I STOP 

IV 

JK 

KNTCT 

L 


Variables, INTERP 

Maximum time (seconds) away from the nearest datum 
that a value VL is interpolated. Zeroes are sub- 
stituted in regions remote from measured data. 
1000} Condensed data from Condenser (see MAIN). 

Time interval (60 seconds) between interpolated 
levels VL. 

DO-loop index, used in initializing VI and T2, 
also used in searching for bracketing COND for 
one-minute interpolation. 

16, 17 "Place markers" which permit starting the search 
for the next bracketing CONDd, ) from the previ- 
ous one. 

1000) Channel and contact number corresponding COND. 
(See DECOM.) 

Pointer (lV-3) for conditional GO TO statement in 
time-bracketing process for interpolation. 

Print file number (File 06) , lOUT ■ 6. 

Terminating condition indicator, ISTOP ■ 2, or » 5 
if contact number LCNTK, or time TBRST is encoun- 
tered, respectively. 

Index indicating variable: contact (pressure), 4; 

Reference, 5; Temperature, 6; or Humidity, 7. 
Length of the filled COND, ICOND array, JK rows. 
The last contact mraber used by INTERP. 

The row index of the one-minute table VL. 
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LCNTK 

LIST 

PCAL( 180) 

Tl(7), T2(7) 

TBRST 

TEST(IO) 

TLPCAL 

TNOH 

Vl(7), V2(7) 
VL(7, 150) 


The last contact nunber calibrated in the PCAL 
table. 

The used length of the VL table. 

The calibrated pressure values at the baroswitch 
contacts. 

Time (seconds from launch) of the bracketing COND, 
ICOND at a given one-minute level L. 

Time (seconds from launch) of balloon burst. 
Diagnostic print control inputs (see MAIN). 
INTER? prints interpolator quantities if TEST(6) > 
O.Ol, e.g., if TEST(6) - 1. 

Time (seconds from launch) that contact LCNTK is 
encountered. It causes termination of data pro- 
cessing. 

Time (seconds from launch) that contact 135 was 
encounterd in DECOM. No Humidity data is pro- 
cessed thereafter. 

Bracketing quantities from COND, ICOND array for a 
given one-minute level L. 

One-minute table produced by INTER? (and by 
ADVANC) . 


VL(1, ) ■ Time (seconds from launch) 

VL(2, ) ■ Azimuth (degrees) 

VL(3, ) ■ Elevation (degrees) 

VL(4, ) ■ Pressure (mb) 

VL(5, ) • Reference frequency (Hz) 
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VL(6, ) ■ Temperature (ordinates) 

VL(7, ) ■ Humidity (ordinates) 
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APPENDIX A 


PROGRAM LIST 
(RAWTNPROC) 

The FORTRAN program list of RAWINPROC given below is included 
for reference. The few handwritten annotations indicate certain 
details assvaned in the main body of the document which, in fact, 
differ from the program version existing and "frozen" at delivery 
of this document. They ere at most minor improvements in that they 
do not affect program performance greatly. Those which may have 
noticeable effect are: 

a. MAIN line No. 55, which extends applicability to negative 
input values of TSTART. 

b. MAIN line No. 155, which corrects a minor error in the 
ultimate values of the temperature and humidity near the 
surface. 

c. ANGLE line No. 39, which would provide smoother angles 
for wind computation. 

d. SEARCH line No. 65, which would more accurately exclude 
short dwells from DECON. 

e. DECOM lines No. 168-169, which eliminate erroneous 
assignment of contact mxnber in certain cases. 

f. DECOM lines No. 326-327, which skip short dwells (as 
intended) . 

g. DECOM line No, 336, which retains, in some cases, the 
Temperature dwell adjacent to a Reference dwell. 

The annotation itt DECOM lines No. 356 and 358, though valid, is not 


b 

b 

neceiaary lince INTERP ignore* repeated contact number*. Other 
annotation* delete *ome of the obeolete (inert) code and update 
aone of the connient*. 
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APPENDIX B 


SAMPLE OUTPUT 
(RAUINPROC and ECC-PRD) 

The primary output of RAvaNPROC is the **input card deck'* 
(File 10, 10 ■ 6) for the concluding Activity No. 3, ECC-PRD. The 
content of the *'one-minute data cards" in this "deck" is listed 
with the label "DECOMMUTATEO OUTPUT AT UNIFORM TIME INTERVALS" 
found near the end of Activity No. 2, below. Detailed output of 
DECOM is listed by time and channel if desired (TCST(7) > 0.01). 
The first eleven ten-point samples of raw data (TIME, FREQ) after 
TSTART are printed to verify proper input. All input card deck 
quantities and computed initializing quantities are routinely 
listed. Certain other output messages indicate status and progress 
of the computation. 

The consequent listed output of ECC-PRD is also included. 
Note that ECC-PRD, used in this "no-ozone" mode, lists zeroes for 
ozone quantities. 

The erroneous Humidity dwell at 23.5 minutes (JX ■ 187) was 
due to the reversed order of source deck cards (DECOM line No. 326, 
327). 

The large Temperature frequency at 24.1 minutes (JK <■ 195) 
appears real. Examination of the raw data (output of HETPASSl) 
will verify the presence of this irregularity in the input data. 
Annotation d). Appendix A, would eliminate this irregular point if, 
in fact, it is composed of a short extreme value followed by a few 
aeconds of noise. The irregularity did not occur when run at the 
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Univ«r«icy using Fils 01 prspsrsd by ths Univsrsicy's countsrpsrt 
to METPA8S1. Ths output shown htrs is ont of thrss tost flights 
run at NASA Usllops Flight Ctntcr on Fsbruary 27, 1981| using 
METPA881 tfhich was then under development. 
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APPENDIX C 


JOB CONTROL DECK 
(RAWINPROC) 


$ 

lOENT 

300700, RAWINPROC 

$ 

OPTION 

FORTRAN 


(object decks— MAIN, ADVANC, ANGLE, TRACK, SEARCH, DECON, INTERP) 

$ 

EXECUTE 


$ 

LIMITS 

05,23K,0,3K 

$ 

FILE 

02,A1S 

s 

FILE 

01,NSTDLB,MLTFIL 

$ 

TAPE 

01,XID,S592 

$ 

ICODE 

IBMF 


(input deck 

for RAWINPROC) 

$ 

IDENT 

300700, ECC-PRD09 


OPTION 

FORTRAN 


(object decks— OZONE, ECC, WODC, OZGRID, RL, TEMPCE, WINDS, ADIR) 

$ 

EXECUTE 


$ 

LIMITS 

05,28K,3K 

$ 

SYSOUT 

08 

$ 

FILE 

05, AIR 

$ 

FILE 

09,A2R 

$ 

TAPE 

07, XIR,,,, SCRATCH 

$ 

ENDJOB 


***E0F 



The job 

control deck above is for the second and third of the 

three-part RAWINPROC system. Documentation for the first part. 

METPASSl, is 

presently under way at NASA Wallops Flight Center (W. 

J. 

Speidel) . 

The middle portion, RAWINPROC, is described in this 


document. The final portion, ECC-FRD, ia operational and document- 
ed as NASA Computer Program 3.0.0700, NASA Wallops Computer Program 
Abstracts. Vol. XXVll. 
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APPENDIX D 


>1LE DESCRIPTION 
(RAWINPROC) 

Pil« 01 — Input tape file containing aucceaaive (0.1-second) 
values, TIME, FREQ, AZ, EL, written by the preceding computer 
routine NETPASSl (see Appendix C) . 

File 02 — Output disc file to the succeeding computer 
routine ECC-PRD (see Appendix C, Appendix A (INTER?), and Input 
Card Deck. 

File 03 — Vestigial file used in program development, not 
used in production runs. 

File 05 — Input file for card input (see Input Card Deck) . 

File 06 — Print file for auxiliary and diagnostic printout 
(Appendix B) . 



